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Abstract 



Inclusive beauty quark production in photon-photon colhsions at CERN LEP2 is con- 
sidered in the framework of the fc^^-factorization approach. Both direct and resolved photon 
^ . contributions are taken into account. The unintegrated gluon distributions in a photon are 

0\ , either obtained from the full CCFM evolution equation or from the Kimber-Martin-Ryskin 

Q [ prescription. The predicted beauty cross section reasonably agrees with the recent experi- 
cn I mental data taken by the ALEPH collaboration. We argue that theoretical and experimental 

studies of the azimuthal correlations in heavy quark production at high energies can serve 
as a crucial probe of the unintegrated gluon densities. 



PACS number(s): 12.38.-t, 13.85.-t 



The problem of beauty quark production at high energies continues to be a subject of 
pointed discussions and intense theoretical studies up to now [1]. First results [2] on the 
6-quark cross section in ep-collisions at HERA were significantly higher than the QCD pre- 
dictions calculated at next-to-leading order (NLO) approximation. Similar observations were 
made in hadron-hadron collisions at the Fermilab Tevatron [3] and also in photon-photon 
interactions at LEP2 [4]. In the latter case the theoretical NLO QCD predictions were below 
the experimental data by three standard deviations. Although the latest measurements [5] 
do not confirm the large excess of the first HERA data over the NLO QCD, the problem is 
not solved so far. The disagrement between the experimental data at the Tevatron and NLO 
QCD predictions was reduced by adopting a special nonperturbative fragmentation function 
of the 6-quark into the i?-meson [60 



^A more exotic solution to this problem was proposed in [7]. 



From our point of view, a more adequate solution was found [8] in the framework of /c^- 
factorization approach [9]. The ^t- factorization approach has also been used for a detailed 
description of numerous experimental data on 6-quark production at HERA [10]. However 
the problem of the 6-quark production in 77 interactions is not solved so far in the /c^- 
factorization approach [11-14]. 

Recently the ALEPH collaboration at LEP2 has presented the result on open beauty 
production cross section in 77 collisions [15]. This is the first published measurement in which 
the lifetime information has been used to identify the heavy flavor in two-photon physical. 
The cross section of the process e+e~ -^ e~^e~bbX has been found to be 5.4 ± 0.8 (stat.) ± 
0.8 (syst.) pb which is fully inconsistent with the previous results quoted by the L3 and OPAL 
collaborations [4], namely 12.8±1.7 (stat.)±2.3 (syst.) pb and 14.2±2.5 (stat.)~5'3 (syst.) pb, 
respectively. In the present note we would like to demonstrate that the ALEPH experimental 
data can be described in the /cT-factorization approach also and to propose an additional 
test to distinguish the different unintegrated gluon distribution functions, which are the main 
ungredient of the /cT-factorization (see, for example, [16]). 

Theoretically, heavy quarks in 77 collisions can be produced via direct and resolved 
production mechanisms. In the direct events, two photons couple directly to a heavy quark 
pair. This contribution is governed by simple QED amplitudes (which are independent of 
the gluon density in the photon). In the resolved events, one photon ("single-resolved") or 
both photons (" double- resolved" ) fluctuate into a hadronic state and a gluon or a quark 
from of this hadronic fluctuation takes part in the hard interaction. At LEP2 conditions 
the heavy quark production via the double resolved processes is highly suppressed [17] and, 
therefore, it will not be taken into account in our analysis. 

The single-resolved contribution to the 77 —>■ bb process is dominated by the gluon 
component of the photon and has the following form in the fc-r-factorization approach: 



dcri_res(77 -^ bbX) 



/^-^^l--^A(x,k|,,^)|A.r(7/ ^ bb)dl.P 



where A.y{x,k'^, fi"^) is the unintegrated gluon distribution in the photon, |A^P(75'* — > bb) is 
the off-shell (i.e. depending on the initial gluon virtuality) matrix element squared, s is the 



total cm. frame energy and a = J ml + pf.exp(y)/yi. The produced beauty quark has 
the transverse momentum p^, rapidity y and azimuthal angle (pb- The initial off-shell gluon 
has a fraction x of the parent photon's longitudinal momentum, the non-zero transverse 
momentum kr (kf. = — fcfn 7^ 0) and azimuthal angle 0. In accord with the fc^-factorization 
prescription [9], the off-shell gluon spin density matrix is taken in the form 

n^y'(^) = ^^- (2) 

In all other respects our calculations follow the standard Feynman rules. The analytic 
expression for the |A^p(7f7* -^ bb) is given in our previous paper [13]. Note that if we 
average Eq. (1) over the azimuthal angle (p and take the limit kf. —>■ 0, we recover the 
well-known formulas corresponding to the leading-order (LO) QCD calculations. 



^The previous measurements by L3 and OPAL collaborations [4] were based on a fiting the transverse 
momentum of leptons with respect to jets. 



The recent experimental data [15] refer to beauty quark production in the e+e" colhsions. 
In order to obtain the corresponding cross sections, the 77 cross sections need to be weighted 
with the photon flux in the electron: 



da{e e -^ e e bbX) = / fy/eixi)dxi / f^/e{x2)dx2da{^'y ^^ bbX), (3) 

where we use the Weizacker- Williams approximation for the photon distribution in the elec- 
tron: 

ttem A + (1 - Xf , gLx , o™2^ / 1 1 



Here aem is the flne structure constant, mg is the electron mass, Q^in = mlx'^/{l — xY and 
QLx = 6GeVMl5]. 

The unintegrated gluon distribution in the photon ^^(x, k^,/i^) can be obtained from 
the analytical or numerical solution of the BFKL or CCFM evolution equations. In order to 
estimate the degree of theoretical uncertainty connected with the choice of unintegrated gluon 
densities, in the numerical calculations we tested two different sets, namely the CCFM [12] 
and KMR [18] ones. First of them was obtained in [12] from the full CCFM equation 
formulated for the photon, and the second one was obtained from the usual (coUinear) 
parton densitieqj using the Kimber-Martin-Ryskin prescription [18]. These distributions 
are widely discussed in the literature (see, for example, [16]). Other essential parameters 
were taken as follows: the 6-quark mass mi, = 4.5 ± 0.1 GeV and the renormalization and 



factorization scale /i = ^Jml + {pr), where (p^) is set to the average pf. of the beauty 
quark and antiquark. In order to investigate the scale dependence of our results we vary 
the scale parameter ^ between 1/2 and 2 about the default value ^ = 1. For completeness, 
we use the LO formula for the coupling constant a;s(/i^) with Uf = 4 active quark flavours 
and Aqcd = 200 MeV, such that agiM"^) = 0.1232. The multidimensional integration has 
been performed by means of the Monte Carlo technique, using the routine VEGAS [20] . The 
full C++ code is available from the authors on requeslo. This code is identical to that used 
in [13, 14]. 

The results of our calculations are displayed in Figs. 1 — 4. Fig. 1 confronts the total cross 
section a{e~^e~ -^ e^e'bbX) calculated as a function of the total cm. energy ^/s with recent 
experimental data [15] taken by the ALEPH collaboration. The solid and dash-dotted lines 
correspond to the results obtained with the CCFM and KMR unintegrated gluon densities, 
respectively. The upper and lower dashed lines correspond to the CCFM gluon density with 
6-quark mass and scale variations as it was described above. Separately shown (as a dotted 
line) is the contribution from the direct production mechanism 77 —>■ bb. It is clear that at 
^/s ~ 200 GeV the cross section is mostly controlled by the single-resolved contribution, i.e. 
jg* -^ bb subprocess. Despite the fact that the central predictions are slightly lower than 
the measured cross section, we observe a reasonable agreement between our calculations and 
the ALEPH experimental data [15] within the theoretical and experimental uncertainties. 
The CCFM-evolved gluon density gives slightly larger cross section compared to the KMR 
one, where the small-x logarithms are not taken into account [18]. A similar effect (but 



■^In the numerical calculations wc have used the standard GRV (LO) parametrizations [19] of the coUinear 
quark and gluon distributions. 
"'lipatov@theory.sinp.msu.ru 



much more clear) has been demonstrated in [10] where, in particular, the beauty photo- and 
lepto-production at HERA has been studied. Note that the sensitivity of our results to the 
variations in the scale n and beauty mass rUb is rather large. However, this sensitivity is of 
the same order approximately as in the massive NLO QCD calculations [21]. 

The transverse momentum and pseudo-rapidity distributions calculated at the averaged 
total e~^e~ energy y/s = 196 GeV (130 < ^i < 209 GeV) are shown in Figs. 2 and 3. As a 
representative example, we have used the following cuts: pr < 20 GeV and \ri\ < 2. In our 
calculations we took into account for both the beauty quarks and anti-quarks. One can see 
again that the difference between the CCFM and KMR predictions is not significant, except 
at large pt (namely pr ~ 10 GeV) only. A similar observation was also made [14] in the case 
of charm production at LEP2. It was shown that the shape and the absolute normalization 
of D* transverse momentum and pseudo-rapidity distributions practically do not depend on 
the unintegrated gluon density. 

We would like to stress that further understanding of the process dynamics may be 
obtained from the angular correlation between the transverse momenta of the produced 
quarks. These quantities are particularly sensitive to high-order corrections. So, in the naive 
LO collinear approximation of QCD, the distribution over A0 = (pb — 4>b must be simply 
a delta function 6{A(f) — tt) since the produced quarks are back-to-back in the transverse 
plane. Large deviations from these values may come from higher-order QCD effects. In 
the /cT^-factorization approach, taking into account the non- vanishing initial gluon transverse 
momentum k-p leads to the violation of this back-to-back kinematics even at leading order. It 
is an illustration to the fact that the LO /cT-factorization formalism incorporates a large part 
of standard (collinear) high-order corrections (see also [9, 16] for more information). The 
differential cross section da/dA(j) calculated at ^/s = 196 GeV is shown in Fig. 4. One can 
see that the shape of this distribution predicted by the CCFM and KMR gluon densities are 
strongly differ from each other. At large A0 ~ tt both gluon densities under consideration 
give similar results, whereas at low A0 ~ the difference is about a factor of 2 in the 
absolute normalization. This fact is directly connected with the properties of non-collinear 
evolution model. Therefore these correlations can be used to constraint the unintegrated 
gluon distributions. A similar effect was also pointed out in the case of beauty production 
at the Tevatron [8]. 

In conclusion, we would like to emphasize additionally that the /c-r-factorization approach 
supplemented with the CCFM-evolved gluon density agrees well with the numerous data on 
the 6-quark production at HERA and Tevatron (without any special assumption on the b- 
quark to 5-meson fragmentation function), as it was demonstrated earlier in [10, 8]. So we 
can conclude that at present there is no contradiction between the CCFM-based theoretical 
predictions and available data on the beauty production at high energies, and we believe 
that the /cT-factorization holds a possible key to understanding the production dynamics at 
high energies. 
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Figure 1: The beauty total cross section a{e^e^ — * e^e^bbX) as a function of the e'^e' 
center-of-mass energy y/s. The sohd and dash-dotted hues correspond to the results obtained 
with the CCFM and KMR unintegrated gluon densities, respectively. The upper and lower 
dashed lines correspond to the CCFM gluon density with variation in 6-quark mass and scale 
as it was described in text. Separately shown is the contribution from the direct production 
mechanism (dotted line). The experimental data are from ALEPH [15]. 
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Figure 2: The differential beauty cross section da/dpr for tlie process e^e 
\ri\ <2 and ^/s = 196 GeV. Notation of curves is tlie same as in Fig. 1. 
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Figure 3: The differential beauty cross section da/d\ri\ for the process e+e -^ e+e hhX at 
Pt < 20 GeV and ^/s = 196 GeV. Notation of curves is the same as in Fig. 1. 
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Figure 4: The differential beauty cross section da/A(f) for tlie process e"'"e~ 
^/s = 196 GeV. Notation of curves is tlie same as in Fig. 1. 
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